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We describe how the fermionic functional renormahzation group (fRG) flow of a Cooper+forward 
scattering problem can be continued into the superconducting state. This allows us to reproduce 
from the fRG flow the fundamental equations of the Eliashberg theory for superconductivity at all 
temperatures including the symmetry-broken phase. We discuss possible extensions of this approach 
like the inclusion of vertex corrections. 
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I. INTRODUCTION 

A considerable part of the current understanding of condensed matter is based on mean-field theories. For super- 
conductivity, the Bardeen-Cooper-Schrieffer (BCS) theory or its counterpart with a realistic phonon spectrum and 
Coulomb repulsion, the Eliashberg theor yii"i i ^ i i i - i, is a powerful framework that allows one to calculate many properties 
of superconductors such as transition temperatures and excitation gaps even quantitativelji^ from a clearly defined 
starting point. Yet in many modern materials these standard pictures seem to be challenged. In each case the question 
arises if the discrepancies arise due to completely new physics or if the established methods fail only partially and can 
be adapted. A mean-field approach can be expected to be good if the fluctuations about it are not too strong and if 
the correlations entailed by the chosen mean-field are far stronger than other types of correlations. Also this second 
requirement is fulfilled less clearly in many low-dimensional or strongly correlated materials, where one often sees a 
competition between various ordering tendencies at low temperatures. Hence, in order to obtain an understanding of 
the situation it is necessary to find a theoretical description which can include all the important channels. Renormal- 
ization group (RG) approaches have the ability to treat a competition between various types of fluctuations in great 
detail. They allow for an unbiased detection and qualitative comparison of Fermi liquid instabilities. This has been 
demonstrated in many works on one-dimensional systemaSiX and the two-dimensional Hubbard niodel^'^'^'^'^^'^^'"'^'^. 

Recently a functional renormalization method has been applied to electrons coupled to phononsi^. In this work the 
authors describe how Eliashberg theory can be reproduced for temperatures at and just above the superconducting 
transition, but not below it. Although this work develops an essentially correct and useful physical picture, it 
highlights a general drawback of the approximate methodsSiSiifliiiiiSiiiii used until now. In these approaches the 
renormalization group flow cannot be extended into the symmetry-broken phase. A notable exception are theories 
where the four-fermion interactions are bosonized at some scalei^, but then one has to work hard to remove the bias 
introduced by that into the flow. In addition, for many problems the type of Hubbard-Stratonovitch decoupling is 
not obvious from the outset. Hence, although several routes seem worth pursuing, a continuation of the flow in the 
unbiased fermionic picture is desirable. In the language of flowing coupling constants in a purely fermionic RG, a 
flow to strong coupling occurs when the interactions seem to diverge at a finite energy scale and consequently the 
perturbative flow has to be stopped. In many cases the flow to strong coupling takes place only in a well defined 
channel of the fermionic interactions and only one specific susceptibility becomes large. This can then be interpreted 
as an indication for spontaneous symmetry breaking in the corresponding channel. The energy or temperature scale 
where this happens can be taken as an upper estimate for the ordering temperature. 

In a previous worlii& we have developed a method that allows us to continue the fermionic RG flow into the symmetry 
broken phase. The idea is to include a small symmetry breaking field in the initial conditions of the flow. This small 
offdiagonal sclfenergy, e.g. a small superconducting gap, grows at the scale where the flow to strong coupling takes 
off but - if one uses a reorganization of the flow equations proposed by Kataninii - prevents a divergence of the 
interactions at nonzero RG scale. This allows us to integrate out all modes down to zero scale and the interactions 
and the offdiagonal sclfenergy saturate at finite values. The Goldstone boson related to the broken symmetry has a 
small mass gap due to the initial symmetry-breaking field. This can be sent to zero afterwards. In Ref. we only 
considered the reduced BCS model with a static attraction and showed how the exact gap value of the BCS theory is 
recovered. Here we generalize the approach to include dynamical phonons and the forward scattering channel. This 
yields an additional equation for the normal sclfenergy. Our extended RG scheme allows us to reproduce the Eliashberg 
equations not only for the symmetric phase above the critical temperature as in Ref. but at all temperatures. 
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Furthermore our treatment makes clear which approximations are used and what one can do to go beyond EHashberg 
theory. 

This paper is organized as follows. In Sectional we briefly introduce the model and the Eliashberg equations. In 
Sec. mil we describe the general fermionic functional RG formalism. In Sec. IIVI we apply an approximate version 
of this formalism to the Eliashberg problem. In Sec. we conclude with a discussion of future extensions of the 
approach. 



II. THE MODEL AND ELIASHBERG THEORY 

The model we study is the basically same as in Ref. We consider spin-1/2 electrons with a dispersion e{k). In 
order to keep the formalism simple we restrict ourselves to a spherically symmetric system with a smooth and finite 
density of states in the energy window of interest around the Fermi level. The electron-electron interaction is given 
by a static and spin-rotationally invariant interaction Vc ~ A:2, fe, A:4). We have written ki for the Matsubara 
frequency and wave- vector of electron i. In our notation, k\ and ki belong to the two incoming electrons and the 
spin component of k\ and the first outgoing particle k-^ is the same. Vc can be envisaged as a screened Coulomb 
interaction. Its precise form is not needed in this paper. Next we add to a phonon-mediated interaction 

^ph(fci,fe,fc3,A:4) = -g{kx,k'i)g{ki,k/^)D{kx - fcs) (1) 

which arises due to the exchange of phonons with propagator 
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Diqo..q) = / dujB{uj,q) 
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.g(fci, k^) is the coupling to the fcrmions. It is well known that the exchange of phonons can induce superconductivity. 
In this case the selfenergy of the electrons acquires a nonzero off-diagonal part A(w) in addition to a quasiparticle 
renormalization factor Z{lu) in the diagonal selfenergy. The Eliashberg theory aims at calculating A{uj) and Z{lli) for 
a given phonon spectrum B{Lo,q). It is contained in a set of two self-consistent equations for the normal self-energy 
and the anomalous self-energy A(fc) which read^ (assuming spin-singlet pairing for simple notation) 

= -Y.VMk,k' ,k' ,k)G{k') (3) 

k' 

A(fc) = -Y^V,s(k,-k,k',~k')F{k') (4) 

k' 

G{k) and F{k) are the diagonal and off-diagonal propagators, respectively, and 

Fcff (fci, /c2, fcs, ki) = u{ki, fc2, fcs, ki) + Fph(fci, /c2, fcs, ki). (5) 

In many cases it is a good approximation to assume that the effective interaction depends only on the frequency- 
momentum transfer ki — k^. The two equations ^ and Q can be visualized in two Fock-type self-energy diagrams 
as shown in Fig. ^ The Z-factor is found from the normal self-energy via 

T.{k)^Y.e{k) + ikQ[l~Z{k)] . (6) 

The even- frequency part Se(fc) is typically weakly fc-dependen1ii. Then it is treated as a renormalization of the 
chemical potential and is not considered further'^ . The spectral gap is determined by the retarded gap function 

In these equations, phonon vertex corrections for the phonon-electron vertex are neglected based on Midgal's theorem. 
The dressing by the Coulomb interaction is absorbed into an effective electron-phonon matrix element. The Coulomb 
vertex is taken as a constant in the static limit^. In this approximation, the Hartree selfenergy contribution is 
ignored. The conventional argument is that the phonon-mediated interaction in this diagram does not produce a 
frequency dependence of the selfenergy and hence no contribution to Z, and possible frequency dependences from 
the renormalized Coulomb repulsion are argued to be less relevant for the questions of interest. However for general 
situations, this diagram can still cause Fermi surface deformations. Hence in problems where the Fermi surface 
geometry and location are important, the Hartree term should be taken into account. 



FIG. 1: Diagrammatic expression for the self-consistent Eliashberg equations. The internal solid lines are full normal (a)) and 
anomalous (b)) propagators. The dashed line denotes the effective interaction. 



III. FUNCTIONAL RG 



The backbone of our approach is the functional rcnormalization group technique for 1-particle irreducible (IPI) 
vertex functionsiSiiS. It is derived from an exact equation for the generating functional for the IPI vertex functions 
of the theory when a flow parameter in the quadratic part Q of the action is changed. Inserting into this equa- 
tion an expansion of the generating functional in monomials of the source fields with the n-point vertex functions 
7*-"^ (iiTi , K2, ■ ■ ■ , Kn) as coefficients one obtains an infinite hierarchy of flow equations for the 71-point vertices. The 
index K = (c, s, /c, iuj) comprises the wave-vector k, Matsubara frequency icj, spin s and the Nambu particle-hole 
index c which is c = + for a fermionic field 'ip[k^iuj, s) and c = — for ijj(~k, —iuj, — s). All bold-face quantities carry 
these Nambu indices. We shall treat the translation invariant case. Then 



(8) 



G{Ki,K2) = S{ciki + C2k2)G{Ki, K2), and similarly for S. Here we used Ki — {ci,si,ki) with fc; — {ki,iuJi). This 
hierarchy of equations is, in a first approximation, truncated after the irreducible 4-point vertex. Then one is left 
with one flow equation for the 2-point vertex, which gives the flow of the self-energy, and one equation for the flow of 
the four-point vertex. They read 



t{K\,K2) = -l 1^^HKi,K2,Pi,P2)S{PuP2) 



(9) 



Pl,P2 



L(Pi,P2,P3,A) 

Pl.P2,P3,Pi 



■ 7(4) {K, , K2, P2,P^W^^ (P4, Pi , i^3, K 

+ 7(4) (A'l, if4, ^2,^3)7^^^ (Pi, Pi, ^2, 

The dot denotes the derivative d/dA with respect to the RG scale which we choose as an infrared cutoff A. S(Ki, K- 
is the so-called single-scale propagatoni^ 



<3) 



(10) 



2) 



S{Ki,K2) = - 5] G{Ki,K)Q{K,K')G{K\K2 



(11) 



with the full scale-dependent propagator G(Ari, ^2). L(Pi, P2, P3, P4) is the scale-derivative of the product of two 
full Greens functions, 

L(Pi,P2,P3,P4) - ^ [G(Pi,P2)G(P3,P4)] (12) 

Equation H12|) corresponds to the modified IPI-RG scheme as introduced by Katanir^ and discussed thoroughly in 
Ref.Ji). This modification is essential in order to obtain correct results in the symmetry-broken phase. The single-scale 
propagator is related to the scale-derivative of the Green's function by 



G{Ki,K2)^S{Ki,K2)+ Y G{Ki,K)t{K,K')G{K',K2) 

K,K' 



(13) 
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FIG. 2: RG equations for the two-point a) and the four-point vertex b). The slashed hne denotes a single-scale propagator Si{p) 
while the dashed line symbolizes the scale derivative of the one- loop diagram. The one-loop graph in b) includes particle-particle 
and particle-hole contributions. With our truncation 7^™'' = for m > 6 the feedback of the 7^''' on 7^*' is neglected. 




FIG. 3: a) Coupling function VA(fci, fcg, ^3). b) Fock diagram for the selfenergy. c) Hartree diagram for the selfenergy. 
d) Particle-particle diagram, e) crossed particle-hole diagram, f) vertex corrections, g) electronic corrections to the phonon 
propagator or screening. 

In the symmetric phase, the spin-rotationally invariant interaction vertex can be expressed^aii by a coupling function 
VA(fci, k2, k^) where the spin indices s (s') of the first (second) incoming and the first (second) outgoing particles are 
the same (see Fig. Oa)). That is, 

7(4) {{+, s, fci), {+, s', fc2), (-, s', fca), (-, s, ki)) = S{ki +k2-k3- ki) VAiki,k2, fca), (14) 

so that VA(fci,fc2,fc3) describes the scattering (fci,s) {k3,s) and (fc2,s') (fc4,s'), and all other values of 7'^^ 
are fixed by the fermionic antisymmetry and the invariance under spin rotations and charge conjugation. With the 
assumptions of the previous section we can write the initial (bare) coupling function as 

VAo(fci,^2,^3) = V;ff(fci, fc2,/c3,^4) = u(fci, A;2, /ca, ^4) - g{ki,kz)g[k2,ki)D{ki - ks) . (15) 

Actually the precise form of the initial interaction is not important, provided it is not long-range, i.e. singular in 
momentum space, and in the following we just study a general regular four-point term Vaq (fci, fe, ^3). In fact, 
the interaction corresponding to an exchange of acoustic phonons can be treated without difficulty even though its 
derivatives in momentum space are unbounded. If one limits the considerations to singlet superconducting pairing, 
spin-rotation symmetry holds at all temperatures. Then the parameterization of the normal interaction vertex can 
still be used in the U(l)-broken phase. In this phase however, anomalous interaction vertices are generated with an 
unequal number of incoming and outgoing outgoing lines. 

Diagrammatically the equations above arc shown in Fig. |21 In this form, the lines in the diagrams do not have a 
direction yet as the fermionic fields still carry the Nambu index. Resolving the Nambu index for the U(l)-symmetric 
normal phase we obtain two diagrams for the self-energy, one Hartree and one Fock term, and 5 diagrams for the flow 
of the interaction. These contain particle-particle and particle-hole diagrams (see Fig. 

IV. COOPER + FORWARD SCATTERING MODEL 

Now we restrict the analysis to a spherical Fermi surface and energy scales that are low compared to the Fermi 
energy. Furthermore we concentrate on systems which do not violate spin-rotation invariance. In the low-energy 
regime the theory simplifies drastically because most contributions are subleading and can be dropped. This allows 
us to reproduce the Eliashberg equations. Let us consider the fiow in the normal phase first. 
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At scale A in the flow, we have already integrated out all degrees of freedom with energies above A. For small 
A, only a thin shell of width A ^ Ep is left to integrate over. The geometry of this thin shell leads to strong 
kinematical constraints if the Fermi surface is convex and positively curved^fiiSiiSSiSi, and there are only three classes 
of two-particle interaction processes where all incoming and outgoing momenta are close to the Fermi surface. These 
are the Cooper processes with zero total incoming wavevector fci + ^2 = 0, the direct forward scattering with fci = fca, 
^2 — and the exchange forward scattering ki — k^, k2 — k^. 

When these momentum configurations are put in as external momenta for the loop contributions in Fig. |21 they 
impose further restrictions on the loop variables as well because of the scale restrictions on the propagators. In some 
of the diagrams, the internal loop variable remains free to explore the full shell around the FS. These diagrams will 
give the leading contributions to the flow in the limit K/ Ep — > 0. They are the a) particle-particle diagrams with 
zero total wavevector with dominant contribution for zero total incoming frequency, therefore we keep VA(fc, —A, k')] 
and h) the particle-hole diagrams with zero wavevector transfer, where the largest contribution (which is possibly the 
only nonzero one) arises for zero frequency transfer. Hence we keep the forward scattering processes VA(fc, A:', k) and 
Vf^ik, k' , k'). The rcnormalization of all other processes not belonging to one of the three classes have less low energy 
phase space and will be suppressed by a factor K/Ep- This can be formalized nicely in a l/iV-expansionSfliSiiSi. 

We emphasize that an additional argument is needed when using this reasoning for the IPI scheme because the full 
propagator, which appears on some of the lines in the diagrams, is not restricted to scale A, but has support on all k 
with |e(fc)| above A. It is the decay of the propagator as one moves away from the Fermi surface that singles out the 
contribution of the above diagrams, where both internal momenta are close to scale A, as leading. 

The steps of reducing the full RG equations to the ones giving Eliashberg theory are as follows. 

• The Hartree terms for the self-energy (diagram c) in Fig. |2l will only give a frequency-independent contribution 
which might reshape the Fermi surface. As we are interested in isotropic systems, the only effect is a constant 
shift which can be absorbed in a redefinition of /x. Thus, for calculating dynamical effects like the Z-factor, we 
may drop the Hartree term. Then the self-energy is only of Fock-type (diagram c) in Fig. O, as in Fig. 

• The particlc-hole terms for the flow of the interactions contain two diagrams which correspond to vertex correc- 
tions (see diagrams f) in Fig. ^ are subleading. These will be dropped. For an interaction vertex corresponding 
to the exchange of acoustic phonons, we can invoke Midgal's theorem, which states that in that case, vertex 
corrections are of order c/vp, where c is the velocity of sound. Note that in this special situation, no small-A 
argument is needed since this holds at all scales. The mixed contributions, those from other phonon types and 
those due to the electron-electron interactions can be argued to be absorbed already in effective vertices^, or, 
more rigorously, be shown to generate only subleading terms in the particle-particle flow by overlapping loop 
estimates^. 

• The bubble particlc-holc diagram (diagram g) in Fig. O corresponds to screening of the Coulomb force or 
to a rcnormalization of the phonon propagator. These effects will already be included in the realistic phonon 
spectrum which is normally used in Eliashberg calculations. Hence we drop these terms as well. 

Hence the remaining diagrams are the Fock term, the particle-particle diagram and the crossed particle-hole ladder 
(diagrams b), d) and e) in Fig|3J). A closer look at the them reveals that the direct forward scattering is renormalized 
by the vertex corrections f) and bubble diagram g) which we have dropped already. Hence, in this approximation, the 
direct forward scattering does not flow in the normal state. This consistent with also neglecting the Hartree diagram 
for the selfenergy, which would contain the direct forward scattering. The exchange forward scattering VA(fc,fc',fc') 
flows with the crossed particle-hole diagram e), which we have kept. Although the contributions in this channel are 
generally small for A T, a rapid flow can developS^ for A ~ T. The exchange forward scattering is the interaction 
occurring in the Fock diagram for the normal selfenergy, which we keep as well. 

A strong flow can also develop in the Cooper channel with zero total incoming wavevector and frequency, possible 
leading to a Cooper instability. The corresponding processes Vi^{k, —k, k') are renormalized by the particle-particle 
diagram c) in Fig. |31 The particle-particle ladder with zero total momentum renormalizes the interaction that appears 
in the anomalous self energy. It is responsible for the growth of the superconducting gap amplitude at the Cooper 
instability. 

Both types of interactions, the Cooper pair scattering (fc, — fc) — > (fc', —k') and the exchange forward scattering 
{k,k') —>■ {k',k), are special in the sense that they are determined by only two wavevectors/frequencies instead of 
three. The exchange forward scattering is described by the vertex functions 'y'"^\K, FK, K' , FK') where the first 

two entries belong to one incoming (e.g., K — {iko,k, s, —)) and one outgoing particle (e.g., FK = (ifco, fc, s', +)). 
FK differs from K by the Nambu index ± and possibly by the spin index s' . For pairing between K and PK {PK 
reverts frequency and wavevector of K , but leaves the Nambu index unchanged, the spin Ps depends on the type 
of the pairing, for singlet pairing it is reversed) the Cooper channel is described by 'y^^\K, PK, K' , PK'). If we 
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FIG. 4: Anomalous interaction vertices generated by anomalous Gorkov propagators and forward+Cooper scattering processes. 
In a) a 4:0 vertex with four incoming lines is generated by two anomalous propagators on the internal lines, in b) one normal 
and one anomalous propagator create a 3:1 vertex with 3 incoming and one outgoing line. The new vertices fall into the same 
category of forward+Cooper scatterings as described by Eq. 1161 



require spin-rotation invariance, the vertices 7^^' {K, FK, K' , FK') and 7^^'' [K, PK, K' , PK'), which do not conserve 
the total spin, are zero. 

Denoting the colhsion partner FK or PK generally by K, we see that the flow due to the two ladder diagrams can 
be expressed in one ladder-type equation for the forward-f Cooper processes. 



^f{K,k,K',k') = ^l j2^(^\K,K,K'\k")-^ [ga(a-, a-")Ga(a-,x")1 j^^\k'\ k'\ K' , K') 



K'' 



(16) 



This means that we drop the two other contributions in Eq. IIUI In those terms K and K would occur each in a 
different vertex of the two 7*^'*^ on the right hand side, hence restricting the low energy phase space on the internal 
lines. 

The nice feature of Eq. H16|l is that it not only describes the flow of the normal 2:2 interaction vertices with 
two incoming and two outgoing lines. It also captures the flow of anomalous vertices which are generated when we 
include an initial gap amplitude AAn(fc) into the flow which breaks the global U(l) invariance explicitly, and which 
is sent to zero after the flow is performed, to induce spontaneous symmetry breaking^^. Furthermore, even if spin- 
rotation invariance is assumed in the normal state, Eq. 1161 can be used to describe situations where the total spin 
is not conserved, e.g. when the superconducting state prefers a certain spin direction. In any case, the new vertices 
generated by Eq. E| still belong to the class of forward and Cooper scatterings and can be described by pairs of 
generalized wavevectors K and K. Therefore Eq. E| remains a closed set and the diagrams dropped above are not 
generated in the flow. For nonzero AAo(fc), the internal propagators also have anomalous (Gorkov) contributions 
which are offdiagonal in Nambu space. Then it is easy to see that ladder diagrams with two anomalous propagators 
(see Fig. 01a)) in Eq. 1161 generate also 4 : or : 4 vertices with four incoming or outgoing lines. The 4 : vertices 
turn out to be essential^^ to stop the flow of the gap amplitude for A — > 0. Moreover diagrams with one normal and 
one anomalous propagator (see Fig. 01b)) lead to 3 : 1 and 1 : 3 vertices with only one incoming or outgoing line. For 
certain situations these vertices can be argued to be absent^ or small, but in general it will be interesting to explore 
the consequences of these anomalous terms. Later we will see that for the special case of the ladder flow, the 4 : 
and 3 : 1 vertices, although nonzero, disappear again out of the final gap equations. This explains why these vertices 
are usually not encountered in standard BCS-type theories. 

The flow of the selfenergy can be expressed with the same class of generalized forward and Cooper vertices. 



SA(i^, k) = -\y. j^j^\k, k, K', k') Sa{k', k') 



(17) 



Here we have assumed that translational invariance on the lattice and in imaginary time is not broken. Therefore only 
normal Nambu-diagonal selfenergies Sii(fcs, ks') and T.22{ks, ks'), and the anomalous Nambu-offdiagonal selfenergies 
T,2i{ks, —ks') or Si2(fcs, — A:s') are nonzero. Eqs. 1161 and [T7I reduce to simple BCS model of Ref. il6| if we set the 
forward scattering to zero and simplify the effective pairing interaction to a constant Vesiki, k^^k^) = —g < 0. 
Thanks to its ladder structure, Eq. ^Ican be solved for any scale A by the Bethe-Salpeter-like equation 



-if{K,k,K',k') = -i'-;^l{KS,K\k' 



(4)/ 



-ffl{K^k,K',k') 



-GA{K'\k")GA{k'\K") 



^f{K'\k",K\k'). 



(18) 



Now we insert this solution into Eq. El This gives 

t^iK^k) = -\Y.-if^{K,k,K',k')^A{K',k') 



K' 
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j^^\k",K",K',K')Sa{K',K') 



+ 1 E l^£iK,K,K",K") 

K',K" 

if' 

= -\j2^^A!iK,K,K',K')GAiK',K'). (19) 

In going from the first to the second equation, we have inserted Eq. 1171 For the last hne, we have used Eq. 1131 and 
translational invariance. We can immediately integrate this equation from Ao down to A = and obtain 

1:aMK,K) - ^Ao{K,K) = -lj2fAliK,K,K',K') Ga=o{K',K') . (20) 

^ K' 

Now we can let 'Ea„{K, K) 0. Resolving the Nambu indices we arrive at two familiar equations. For the normal 
selfenergy I](fc,s) — 'Eii^a=o{K, K) for K — k,iko, s, + and K = fc,ifco,s',— we get, picking the Fock contribution 
(and again assuming spin-singlet pairing for simplicity), 

m) = -T.^Ao{k,k\k')G{k'). (21) 
fc' 

The factor 1/2 goes away in the summation over the internal Nambu indices. Defining the offdiagonal propagator 
F{k) — G2i(fc, s — fc, —s) and the gap function A(fc) — I]2i(fc, s, — fc, — s) we also get the second selfconsistency equation 

A{k) = -Y,VAo{k,-k,k')F{k'). (22) 
fc' 

These equations are equivalent to Eqs. |31and0] Therefore the Eliashberg theory is contained in the RG approach when 
the flow is restricted to the ladder-type diagrams in Cooper and exchange forward scattering. The anomalous vertices 
with unequal numbers of incoming and outgoing lines have disappeared again from the equations. For establishing 
the connection to the selfconsistent equation we did not have to calculate them at all. The equivalence of ladder-type 
RG equations to selfconsistent equations was shown for the normal state by Katanin^i, who pointed out that in this 
case the RG flow fulfills the respective Ward identities, and for the superconducting state in Ref.^^ In our case the 
Ward identity for the global U(l) invariance holds and the emerging Goldstone boson can be identified in the flow of 
the interactionsiS.. 

Note that after confining the analysis to the separate channel for forward and Cooper scattering, respectively, no 
further assumptions have gone into the derivation. In particular, no simplification of the frequency or wavevector 
dependences are necessary. Furthermore, the approximation of the normal selfenergy in terms of a Z-factor is not 
essential to establish the correspondence of the renormalization group and the self-consistent Eliashberg formalism. 
We also note that the ladder equation (|16|l and the selfenergy equation Ijl?!) already go beyond Eliashberg theory. 
Eq. El also contains the bubble summation for the direct forward scattering (fcs, fc's') {ks,k's') of the diagrams 
g) in Fig. 13 and the corresponding Hartree selfenergy, which is diagram c) in Fig. |31 These terms are usually not 
considered in Eliashberg theory, for the reasons indicated above. We have left out these terms in the RG as well, in 
order to exhibit the approximations involved in obtaining Eliashberg theory most clearly. Including them is analogous 
to the procedure for the exchange forward scattering and easy because it just adds a Hartree term to Eq. [TTI 

Let us briefly discuss the wavevector and frequency dependences and how they develop during the flow from Aq to 
A = 0. For this we assume that the initial pairing interaction is factorizable in the incoming and outgoing indices, 

^Aoik, k') = E m{k)gdk') (23) 

with orthonormal functions gt{k) in wavevector-frequency space. Then it is not difficult to see that with the assumed 
ladder structure of the flow in the Cooper channel, the RG equations decouple into separate equations for the coeffi- 
cients Vi^A- Positive Vi^a will decrease during the flow. This also includes the s-wave repulsion due to the Coulomb 
interaction as described by Morel and AndersonS^. It will be interesting to see if additional insights about this process 
can be learned with the RG approach. Negative Vg^A will increase until a Cooper instability is reached the most 
negative ^-channel at a critical scale A^. At this point the normal pairing interaction VA{k, —k, k') becomes large. 
Above this scale, if we start with a small initial gap amplitude Aq, i.e. small off-diagonal propagators, the additional 
anomalous couplings and the gap amplitude AA(fc) will have departed only little from their initial values. However 
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at they develop a strong flow due to the almost-divergence of VA(fc, k'). Now it is clear from Eqs. 1161 and [T7I that 
the fc-dependence of the anomalous couplings and A/^(k) follows the /c-dependence of the divergent component in the 
Z-channel, given by gi. In fact, if we start with an initial gap Ao(A:) = J2e ^o,e9e{k), only the /-component with a 
Cooper instability at this scale will be pulled up and converge to a nonzero value for A ^ in the limit Ao(fc) 0. 
All other ^-components in the final gap amplitude Aa=o(A:) without a Cooper instability in their channel disappear 
for Ao(fc) 0. The divergence of the interactions is stopped just below A^ by the rapid growth of the gap. The 
final value of the normal interaction vertex scales cx I/Aq.^ . In practice this bound can be used to keep the couplings 
near the perturbative range. The convergence of the results in the limit Ao(A:) — > has been discussed in Ref. Il6l 
All other ^-components in the final gap amplitude Aa=o(^) without a Cooper instability in their channel disappear 
for Ao(fc) — > 0. The resulting gap function corresponds exactly to the one obtained by solving the Eliashberg gap 
equation 1)22(1 with the most attractive component Vi only. It will be interesting to see how subsequent transitions 
due to subdominant pairing channels are described in this scheme. 



V. DISCUSSION 



We have shown that a suitable truncation of the fermionic functional renormalization group serves as a basis for the 
Eliashberg equations both in the normal and symmetry-broken phase, thus extending the result of Ref. 14 to temper- 
atures below the critical temperature. In our view our rederivation of the Eliashberg equations from an approximation 
to the RG is of methodical interest — we have described how, by dropping terms, the RG equations reduce to a system 
that integrates to give the self-consistency equations. Namely, integrating the flow of the interactions and selfenergies 
in the ladder approximation described above from the initial bandwidth Aq down to A = produces the same answer 
as solving the selfconsistent equations Q and which the Eliashberg theory is based on. For many cases the latter 
procedure will still be easier because one can work with bare (effective) interactions and does not have to keep track 
of running coupling functions. The valuable advantage of the fermionic functional RG approach presented here and 
in Ref. [Ij is, however, that it paves the way for the inclusion of a number of effects beyond Eliashberg theory. This 
has not been done yet, but while more work is underway let us briefly outline some straightforward extensions. 

Note that the Eliashberg equations (|21|l and H22I) could be derived only under the assumption of ladder-like effective 
vertices in the forward and Cooper scattering channel. In contrast with that, the flow equations ^ and (|10|l hold more 
generally. The difference becomes particularly clear when we consider situations where the initial interaction does 
not contain any attractive component in the pairing channel. Then Eliashberg theory and the corresponding ladder 
approximation to the flow will not flnd a superconducting solution. However, it is known since the work of Kohn 
and Luttinger^^ that particle-hole corrections to the pair scattering will always create an attractive component in the 
effective interaction. In our case, in Eq. ^1 the only coupling of the particle-hole processes into the pair scattering 
is due to the overlap between forward and Cooper scattering which gets small due to phase space restrictions. This 
last argument applies only at very low scales A, and it is in neglecting the attractive effective interaction generated 
by the integration of flelds at higher scales that these approximations miss out the Kohn-Luttinger effect. The more 
careful analysis of the RG flow, which takes the higher scales (regime 1 in Ref. [l^ into account correctly, includes 
the perturbative corrections to the pair scattering, in principle to all orders in the initial interaction, hence renders 
a correct picture of the existence and relative strength of pairing instabilities that are dynamically generated. For 
example, in the two-dimensional Hubbard model on the square lattice near half-flUing it is known that particle-hole 
corrections generate an attractive interaction in the c?2.2„j,2-wave channel^^, and this becomes especially clear in RG 
studies2i2ii&. The extension of the flow beyond the Cooper instability opens a way to analyze the gap structure of the 
d-wave pairing in more detail. 

In addition to the detection of pairing channels that are not present in the initial interaction, the full one-loop flow 
allows one to study the influence of vertex corrections due to the electron-electron interactions on the pairing in a 
systematic way. For example, the particle-hole channel renormalizes the s-wave pairing channel as well. This effect, 
basically due to the same diagrams as the Kohn-Luttinger effect, is known to have quantitative consequences even for 
small interactions^^. Furthermore Migdal's theorem, although proven for acoustic phonons, has been argued to break 
down under various other circumstances^S. 

We emphasize that the concepts described above are not limited to symmetry breaking in the superconducting 
channel. Other types of flows into long-range ordered states can be performed as well. In the context of low- 
dimensional electron systems, the analysis of the interference between magnetic and superconducting order seems 
another promising route. 

Finally we address the interesting issue which types of runaway flows can be brought to a safe end by the inclusion 
of some symmetry breaking field in the initial conditions; this also indicates some limitations of this method. By 
adding a static offdiagonal part to the selfenergy one assumes that that there is long-range order in the system, 
and this is a rather strong assumption. In an exact treatment, the behaviour of correlations at various scales has to 
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be described by a time- and space-dependent order parameter field, while the above ansatz corresponds to a field 
independent of space and time. In low-dimensional systems, and under very general conditions, long range order 
gets destroyed by long-wavelength fluctuations of the order parameter fields, so that long range order gets replaced 
by a Kosterlitz-Thouless like phase for two-dimensional superconductors or phases with even faster decay of the 
correlations of the order parameter field, if they are described by nonlinear sigma models at low energies. In certain 
one-dimcnsional models like Hubbard ladder systems^, the RG often flows to strong coupling but all correlations 
remain short ranged even for T —^ 0. The absence of long-range correlations in this case is due to a combination of 
the above-mentioned effects and the competition of different interactions. One may expect the above ansatz of a static 
symmetry-breaking component to extend to a regime where the dynamics of the order parameter fields takes place 
on much larger temporal and spatial scales than that of the fermions, so that one can hope to describe some aspects 
of phases without long range order. For the understanding of the vicinity of the transition, and the situation where 
strong fluctuations persist to very low scales, one needs to generalize this scheme by including the order parameter 
fields themselves in the description. 
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